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ABSTRACT 

Context. High spectral resolution X-ray observations of classical T Tauri stars (CTTSs) demonstrate the presence of plasma at T ~ 
2-3x10* K and ~ 10" - 10'^ cm"^, unobserved in non-accreting stars. Stationary models suggest that this emission is due to shock- 
heated accreting material, but they do not allow to analyze the stability of such material and its position in the stellar atmosphere. 
Aims. We investigate the dynamics and the stability of shock-heated accreting material in classical T Tauri stars and the role of the 
stellar chromosphere in determining the position and the thickness of the shocked region. 

Methods. We perform 1-D hydrodynamic simulations of the impact of the accretion flow onto chromosphere of a CTTS, including 
the effects of gravity, radiative losses from optically thin plasma, thermal conduction and a well tested detailed model of the stellar 
chromosphere. Here we present the results of a simulation based on the parameters of the CTTS MP Mus. 

Results. We find that the accretion shock generates an hot slab of material above the chromosphere with a maximum thickness of 
1.8x10' cm, density ~ 10" - 10'- cm"^ , temperature T ~ 3 x 10* K and uniform pressure equal to the ram pressure of the accretion 
flow (~ 450 dyn cm"^). The base of the shocked region penetrates the chromosphere and stays where the ram pressure is equal to 
the thermal pressure. The system evolves with quasi-periodic instabilities of the material in the slab leading to cyclic disappearance 
and re-formation of the slab. For an accretion rate of ~ IO^'^Mq yr"', the shocked region emits a time-averaged X-ray luminosity 
Lx « V X 10^' erg s"', which is comparable to the X-ray luminosity observed in CTTSs of the same mass. Furthermore, the X-ray 
spectrum synthesized from the simulation matches in detail all the main features of the O vni and O vii lines of the star MP Mus. 
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1. Introduction 

Pre-main sequence stars are strong X-ray emitters, with X-ray 
luminosity up to lO"' times the solar value (Lx ~ 10^^ erg s"'). 
Similarly to main-sequence stars, the X-rays are likely produced 
by low density (rig ~ 10'" cm"-*) plasma enclosed in coro- 
nal loop structures and heated at temperature T ~ 10^-10^ 
K jFeipelson &Montmerlj|l999l) . During the last five years, 
high resolution (R ~ 600) X-ray observations of some classi- 
cal T Tauri stars (CTTSs) (TW Hya, BP Tau, V4046 Sgr, MP 
Mus and RU Lupi) have shown the presence of X -ray plasma at 
T ~ 2 - 3 X lO'' K and denser than 10" cm"^ dKa stner et alj 
20021: ISchmittetalJl2b0 5; Gun ther et alj 120061: [Argiroffi et al. 
200 7*: Robra de & Schmittil2007i) . which suggests an origin dif- 

fer ent from the coronal o ne. 

ICalvet & Gullbrind ( Il998l) and iLamzinl ( Il998l) proposed 
that. X-ray emission from CTTSs could also be produced by the 
accreting material. In fact the material heated by the accretion 
shock at the base of the accretion column could reach a temper- 
ature T ~ 10* K. Starting from the hypothesis of the presence 
of a stationary strong accretion shock an d solving the c onser- 
vation laws for the shock-heated region, iGtinther et alj (l2007h 
have quantitatively demonstrated that the accretion shock model 
well explains some non coronal features of the X-ray observa- 
tions of TW Hya. However, these models are based on several 



approximations and do not include a detailed model of the stel- 
lar atmosphere, which indeed determines the shock position and 
could influence the profile of pressure and density and, there- 
fore, the thickness of the shocked region. They assume stationary 
conditions, but different studies proved the existence of thermal 
instabilities in presence of a n acc reting flow impacting o nto a 
stellar surface (Laiig er et alJll98U IChevalier & Imamural [19821: 
Koldoba et al. 20081) . Considering that the physical structure of 
the shock-heated region and of its temporal evolution determines 
the fraction of the accreting energy re-emitted in the X-ray band 
and its spectral behavior, further investigations on the issues dis- 
cussed above are required to understand better the X-ray emis- 
sion from shock-heated plasma, to check their observability and, 
then, to derive the physical properties of the accretion process 
from high resolution X-ray observations. 



Here we address these issues with the aid of a time- 
dependent hydrodynamic numerical model describing the im- 
pact of an accretion stream onto the chromosphere of a CTTS. 
As a first applic ation, we compared the X-ray spectra of the 
CTTS MP Mus dArgiroffi etalj 12007 ) with those synfliesized 
from the results of a simulation tuned on this star 
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2. The model 

We assume that the accretion occurs along a magnetic flux 
tube linking the circumstellar disk to the star We consider 
the flux tube analogous to closed coronal loops observed on 
the Sun; therefore, as done for standard c oronal loop mod- 
els (e.g. iPeres et al.lll982t iBettaet al.lll997h . we assume that 
the plasma moves and transports energy exclusively along the 
magnetic field lines. This hypothesis is supported by the low 
value of the plasma parameter p - P/iB^/^n) ~ 10"^, where 
P ~ 4 X 10^ dyn/cm^ is the expected post-shock zone pressure 
and B ~ 10^ G is the typ ical magnetic field at the surface of a 
CTTS (lJohns-Krullll200'i) . Note that plasma density and veloc- 
ity are expected to vary across the section of an accretion stream 
(e.g lRoma n ovaet al . 2003). However, in the case of /? <k 1, the 
stream can be considered as a bundle of "elementary" streams, 
each characterized by different values of density and velocity. 
Our 1-D model describes one of these elementary streams. We 
limit our analysis to the impact of the accretion stream onto the 
chromosphere and consider the portion of the flux tube close to 
the star. We assume the magnetic field lines perpendicular to the 
stellar surface and a plane parallel geometry. 

The impact of the accretion stream onto the chromosphere is 
modeled by numerically solving the time-dependent fluid equa- 
tions of mass, momentum, and energy conservation, taking into 
account the gravity stratification, the thermal conduction (in- 
cluding the effects of heat flux saturation) and the radiative losses 
from an optically thin plasma: 



5p 5pv ^ 
dt ds 



dpv d{P + pv^) 



dpE d(pE + P)v dq 

-T- + z = P^S + Eh- -r-- neriHMT) , 

at as OS 
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where t is the time, s the coordinate along the magnetic field 
lines, p the mass density, v the plasma velocity, P the thermal 
pressure, g(s) the gravity, T the plasma temperature, E - 
the total gas energy per unit mass, e the inte rnal ene rgy per unit 
mass, q the heat flux in the formulation of 'Spitzer ( 1962) cor- 
rected for the effect of flux saturation (Cowie & McKee 1977), 
He and tin are the electron number density and hydrogen num- 
ber density, respectively, A{T) is the radiative losses per unit 
emission measure, f5(p, T) the fractional ioniz ation jn^ln-n) de- 
rived from the modified Saha equation (see lBrown|[l973l for 
the solar chromosphere condition), y- 5/3 is the ratio be- 
tween the specific heats, - 2.14 x lO"^"* g is the average 
atomic ma ss (assuming heavy eleme nts abundances 0.5 the so- 
lar values: lAnders & Grevesselll989h . and E^ is a parametrized 
chromospheric heating function (Eh = for T > 8 x 10"^ 
K) defined, as in Peres et al. (1982), to keep the unperturbed 
chromosphere in stable equilibrium. The radiative losses per 
unit emission measure are derived by the PINTofALE spec- 
tral code (Kashvap & Drake 2000) with the APED VI. 3 atomic 
line database (iSmith et al.ii2001i) . and assuming heavy elements 
abundances 0.5 the solar values ( Anders & Grevesse 19 89), as 
obtain ed from X-ray observations of CTTSs by iTelleschi et al I 
(l2007h . The equations are solved numerically using the flash 



code (iFryxell et al.ll2000l) . an adaptive mesh refinement multi- 
physics code for astrophysical plasmas. The code has been ex- 
tended with additional computational modules to handle the ra- 
diative losses, the thermal conduction, and the evolution of frac- 
tional ionization (ne/ne)- 

The computational domain extends over a range D = 4.34 x 
10^ cm above the stellar surface. We allow for 5 levels of re- 
fine ment in the adap tive mesh algorithm of flash (paramesh; 
Ma cNeice et al. I l2000l) . with resolution increasing twice at each 
refinement level: As^ax =2.1x10^ cm at the coarsest resolution, 
and As,„i„ = 1.3 x 10^ cm at the finest level, which would cor- 
respond to a uniform mesh of ~ 30000 grid points. We analyzed 
the effect of spatial resolution on our results by considering two 
additional simulations which use a setup identical to the one dis- 
cussed here, but with either 4 or 6 levels of refinement. We found 
that the adopted resolution is the best compromise between ac- 
curancy and computational cost and that the system evolution is 
well described in its detail. 

The simulation presented here covers a time interval of about 
2000 s. We used the accretion parameters (velocity and den- 
sity) derived by Argiroffi et al. (2007) to match the soft X-ray 
emission of MP Mus. We calculated the gravity considering the 
star mass M = l.lM p) and the star radius R - 1 .3Rq used by 
lArgiroffi et al.l(l2007l) . 

The external part of the initial configuration, extending from 
2.75 X 10^ to 4 X 10^ cm, consists of an accretion stream con- 
stant in density {rie = 10" cm"^), temperature (T = lO-' K) 
and velocity (v = 450 km/s). The inner part of the initial con- 
figuration consists of a static chromosphere. We reproduce the 
pressure gradient of a young stellar chromosphere by consider- 
ing the temperature profiles prescribed by Vernazza et al.l (Il973h 
solar chromosphere models scaled in order to match a pressure 
of ~ 7 X 10"^ dyn cm"^ at the base of the chromosphere. As 
boundary conditions we consider fixed values both at the top 
(n^^^ = 10" cm-^ Tacc = 10^ K, Vacc = 450 km/s) and at the 
base (nchi- = 10'^ cm-\ T^ur = 4.44 x 10^ K, Vch,- = 0.0 km/s) 
of the computational domain. In principle these boundary con- 
ditions lead to accumulation of matter at the base of the chro- 
mosphere. However, we have estimated that this effect becomes 
significant on a timescale 200 times longer than that explored 
by our simulation. In fact, we checked that, for i < 10*^ cm, the 
chromosphere remains virtually unperturbed (with variations of 
mass density below 1%) during the timescale considered. Note 
also that we neglect the heating of the chromosphere (in particu- 
lar, at the lower boundary) due to the X-ray emission originating 
from hot plasma. In the case of MP Mus (effective temperature 
Teff ~ 5000 K), this approximation is justified by the low ra- 
tio between the energy flux coming from the accretion flow and 

from the photospheric emission, (pa(;ci'acc)/(4o"7^e(f) ~ 0.1, where 
cr is the Stefan-Boltzmann constant. The stability of the chromo- 
sphere was tested by dedicated simulations longer than 100 ks, 
some of which considering also strong transient heating. 

3. Results 

The impact of the accretion stream onto the stellar chromosphere 
leads to the formation of transmitted (into the chromosphere) 
and reverse (into the accretion column) shocks. The latter prop- 
agates through the accretion co lumn building up an hot slab. As 
pointed out by [Arg iroffi et al] (|2007|) . in the strong shock limit 
(IZel'Dovich & Raizen 1 1967b . the expected temperature of the 
slab is ~ (3/16)(//mH/^b)Vflcc ~ 3 x 10^ K and, since it is sub- 
ject to the radiative cooling, the expected maximum thickness is 



where v. 



ps 



[;/4 « 100 km s is the post-shock 
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Fig. 1. Evolution of plasma temperature (A), density (B), pres- 
sure (C), and velocity (D) distributions along the flux tube from 
the chromosphere to the unperturbed accretion stream, sampled 
every 60 s from to 300 s (left panels) and every 30 s from 330 
to 420 s (right panels). The figure shows the inner portion of the 
spatial domain, including the chromosphere and the hot slab. 



plasma velocity in the slab and Trad is the cooling time for the 
shocked gas defined as 



1 P T^^^ 

Trad = T T"^ ~ 6.7 X 10^ 



(4) 



where, for temperatures characteristic of the slab {T w 10^ - 
10^ K), we approximated the cooling function as MT) x 6.0 x 
JQ-20J--I/2 gj-g g-i ^jj^3 Pqj. j-jjg typical values of T and «e of 

the hot slab (7 « 3 x 10*' K and «e = 4«acc ~ 4 x 10" cm^^), 
Trad ~ 90 s and the thickness of the slab is expected to be Dsi^b ~ 
10^ cm. Note that deviations from equilibrium of ionization may 
be present beyon d the shock su rface due to the rapid growth of 
temperature (e.g. iLamzinlll998l) . However, these deviations are 
expected within a distance (/nei ~ VpsTc^y ~ 10^ cm from the 
shock, where Teg, the ionization equilibrium timescale, has been 
estimated for the characteristic density and temperature of the 
slab. Being c/nei ^ i'siab, the non-equilibrium ionization effects 
are not relevant for the evolution of the hot slab. 

The evolution of temperature, density, pressure and velocity 
during the first part of our simulation is shown in Fig. [1] During 
the first 300 s of evolution, the reverse shock gradually heats up 
the accretion column at T ~ 3 x 10^ K and the inflow piles up 
material creating a hot slab with density 11^.-3 



= 450 dyn cm ^ (see left panels of Fig.[T]). The 



base of the shocked region initially penetrates the chromosphere 
and comes at rest where the ram pressure is equal to the thermal 
pressure. The ram pressure, therefore, determines the position of 
the slab inside the chromosphere, while the chromosphere below 
the shock region remains virtually unperturbed. Due to pile up of 
material at the base of the shocked column, the radiative losses 
gradually increase there. At the end of this phase, the extension 
of the hot slab is « 1.8 x lO** cm. 

In the subsequent phase of evolution, the strong radiative 
cooling of the high density plasma accumulated at the base of the 
hot slab triggers thermal instabilities there (Field 1965 ; see right 
panels of Fig. [1): plasma temperature and pressure decrease by 
more than two orders of magnitude in few seconds, leading to the 
collapse of the upper layers of the accretion column. As a con- 
sequence, the residual hot slab fells down and the reverse shock 
moves downward to the chromosphere. The additional compres- 
sion due to the collapse leads to the further increase of plasma 
density and radiative cooling at the base of the hot slab which 
gradually cools down in ~ 100 s. 

After the accretion column is completely cooled, a new hot 
slab is generated (see the dashed-dotted lines in the right panels 
of Fig. [U and the system starts a quasi-periodic evolution with 
alternating heating and cooling phases lasting ~ 400 s linked to 
the up and down displacement of the reverse shock. Our simula- 
tion tracks the system evolution for 5 periods without significant 
differences among them. Therefore, we presume that the results 
can be extended to longer time intervals. Note that the period 
of the fluctuation depends on the parameters of the inflow and 
partially of the chromosphere. The assumption of steady flow 
is appropriate over the limited time explored. The whole issue 
should be revisited for simulations covering timescales compa- 
rable to those of changes due to stellar rotation, magnetic field 
evolution, etc. . 

Note that iLanger et al.l (Il981h have found periodic variations 
of the accretion shock position and, therefore, of the thickness 
of the hot slab in the context of compact objects. As in our case, 
the cycle consists of time intervals when the shock moves up- 
ward and the gas accumulates behind it and phases during which 
the shock moves downward as the hot gas c ools radiatively. In 
the context of CTTSs, Kold oba et al.l (l2008h found oscillations 
of the accretion shock position with periods ~ 0.02 - 0.2, due 
to the same mechanism described here. These short periods are 
mainly due to the set of values of accretion flow parameters (ve- 
locity and density) and heavy elements abundances adopted by 
these authors. In addition, our model includes the eff'ects of ther- 
mal conduction (important in the energy budget) and the stellar 
chromosphere (which determines the shock position). 

We checked the observability of the X-ray emission pro- 
duced by the shock-heated material by synthesizing the spec- 
trum in the energy range [0.5 - 8.0] KeV from the simulation 
results, and by calculating the time-average luminosity over the 
time interval covered by the simulatiorQ Since our model is 1- 
D, we have to assume the flux tube cross section. Considering 
the velocity of the accretion stream v - 450 km s"\ its den- 
sity «e = 10" cm"-' and an accretion rate ~ 10""* Mq yr"\ 
the accretion stream cross section is ~ 6.5 x 10^" cm^. Using 
this cross section, we obtain an X-ray luminosity varying from 
~ 8.4 X 10^' to ~ 1 .6 X 10^"' erg s"' with a time-average value of 



' High resolution X-ray spectra of CTTSs are gathered with exposure 
times of ~ 10 - 100 ks, i.e. much larger than the time interval covered 
by our simulation. 
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Fig. 2. Observed X-ray spectrum of the star MP Mus (gray line) 
with the synthetic spectrum derived from the simulation (black 
line). 



- the thickness of the slab grows until it reaches ~ 1.8x10^ cm, 
when it undergoes a thermal collapse due to the enhancement 
of density and, therefore, of radiative losses at its base; 

- the system oscillates quasi-periodically between the forma- 
tion of an hot slab and its thermal collapse by thermal insta- 
bility; 

- according to our model, an accretion rate of ~ lO^'^Mgyr"' 
is enough to produce an X-ray luminosity varying between 
~ 8.4x10^' and 1.6x10^" erg with a time-average value 
of Lx ~ 7 X 10^^ erg s"', which is comparable to the overall 
X-ray luminosity observed in the CTTSs of the same mass; 

- our model is able to explain in detail the observed X-ray 
spectrum of the star MP Mus ( Argiroffi et al. 200 7). 

A large set of simulations, exploring the domain of the phys- 
ical parameters of the system, will be studied in a future paper. 



Lx = 7.2 X 10^^ erg s"'. Since this value is comparable with the 
typical overall l uminosity observed in the young T Tauri stars of 
the same mass dPreibisch et al.ll2005l) . our model shows that the 
shock-heated material can contribute to the X-ray emission of 
the CTTSs. 

In Fig. |2l we compare the X-ray spectrum synthesized from 
the simulation (averaged over the time interval covered by the 
simulation) with the spectrum of the star MP Mus observed with 
the Reflection Grating Spectrometers (RGS) on board the XMM- 
Newton satellite (Argiroffi et al. 2007). The synthesis of the X- 
ray spectrum takes into account the instrumental response of 



XMM-RGS and the interstellar absorption (A^h = 5x10^" cm"-^) 
derived from the observations, and assumes the emitting source 
at the same distance of MP Mus (86 pc). We also assumed that 
the time-average X-ray luminosity produced by the shock-heated 
slab in the 18 - 23 A wavelength range is equal to the observed 
MP Mus luminosity in the same range. With the above assump- 
tions, it turns out that the accretion rate is ~ 7.7 x 10"" Mq 
yr"' and the accretion stream cross section is ~ 5 x lO^'* cm^ 
(i.e. a surface filling factor of about 0.5%). These values are in 
agreement with t hose (~ 5 x 10~" M fr, yr"' and ~ 3 x 10^° 
cm^) derived by lArgiroffi et al.l (l2007l) from the analysis of the 
observations. 

Figure |2] shows O viii Lya and the O vii triplet. The ratio 
between the O vm Lya and the O vii resonance line is a tracer 
of plasma temperature and the ratio between the O vii forbid- 
den and intercombination lines is a tracer of plasma density. The 
good agreement between the observed and the synthetic line pro- 
files demonstrates that our model is able to explain the origin of 
the whole dense X-ray plasma observed in CTTSs. Furthermore, 
our hypothesis on heavy elements abundances is supported by 
the agreement of the ratio between lines and continuum inten- 
sity. 



4. Conclusions 

We model the impact of an accretion flow on the chromosphere 
of a young T Tauri star. Our main results are: 

- the impact is able to build up a slab of accreting plasma with 
density ng ~ 5 x lO" cm"^ and temperature T ~ 3 x 10^ K; 

- the pressure along the slab is uniform and equal to the ram 
pressure of the accretion shock; therefore, the slab is rooted 
down in the chromosphere where the thermal pressure is 
equal to the ram pressure; 
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